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Abstract 

A comprehensive model study of Bose-Einstein correlation radii in heavy ion 
collisions is presented. The starting point is a longitudinally and transversally 
expanding fireball, represented at freeze-out by an azimuthally symmetric emission 
function. The freeze-out temperature is allowed to feature transverse and temporal 
gradients. Their effects on the correlation radii are studied. In particular, we 
evaluate numerically their dependence on the transverse mass of the particle pairs 
and check a recent suggestion, based on analytical approximations, that for certain 
reasonable source parameters all three correlation radii satisfy simultaneously a 
I/VMT scaling. 



1 Introduction 



Recently, two-particle intensity interferometry, exploiting the effects of Bose-Einstein 
(BE) symmetrization on the two-particle momentum spectra, has been developed into a 
powerful tool for measuring not only the space-time dimensions of the particle emitting 
object, but also its dynamical state at particle freeze-out. Modern Bose-Einstein inter- 
ferometry thus goes far beyond the original work by Hanbury Brown and Twiss, who 
introduced and successfully demonstrated photon intensity interferometry for static 
sources in astrophysics Q, and the pioneering paper by Goldhaber, Goldhaber, Lee, 
and Pais Q who first exploited similar ideas in particle physics. Good reviews of the 
basic theoretical and experimental techniques can be found in ||3|, |5|, |6|, ^ while the 
more recent theoretical developments are summarized in a series of lectures published 
ini. 

Collective dynamics of the source leads to a characteristic dependence of the two- 
particle correlation function on the total momentum of the particle pair. The width of 
the correlator as a function of the relative momentum of the two particles measures the 
size of the "regions of homogeneity" inside the source across which the momentum 
distribution changes sufficiently little to guarantee measurable rates for pairs with 
similar momenta of the two particles. For rapidly expanding sources widely separated 
points do not emit particles of similar momenta; the homogeneity regions seen by 
BE correlation function thus form only a fraction of the whole source. Moreover, the 
momentum of a particle defines the homogeneity region from which it is coming; we 
will call the homogeneity region which contributes to the emission of particles with a 
given momentum the "effective source" for such particles. 

The size parameters of the "effective source" (which are measured by the correlator 
and which we will thus call "correlation radii" or simply "HBT radii" ) are related to the 
magnitude of the velocity gradients in the source, multiplied with a thermal smearing 
factor jM^ originating from the local momentum distribution jl^, |l2|. They 
are thus affected by both the expansion velocity profile of the source and by variations 
of the width ("temperature") of the local momentum distributions across the source 



|11, 13, 14, Hq, M. Another physical mechanism which can, even in the absence of flow 



and temperature gradients, lead to dramatic differences between the "effective source" 
of particles with a fixed momentum K and the momentum-integrated distribution of 
source points in space-time is strong reabsorption or rescattering of the particles in the 
interior of the source; this leads to a surface-dominated distribution of last interaction 
points with a strong preference for outward-directed momenta. Such "opaque sources" 
and their effects on the correlation radii have recently been discussed in |17, 18| on 
a semianalytic level using simple models. Our investigations on this subject will be 
published in the next paper |p^ . 

While the basic relations between such source features as mentioned above and the 
properties of the measurable two-particle correlation functions appear to be qualita- 
tively understood, this is not sufficient for a quantitative interpretation of two-particle 
correlation data and a complete reconstruction of the source from the measurements. 
Quantitative numerical investigations, including comprehensive parameter studies, so 
far exist mostly for transparent expanding sources with constant freeze-out temper- 



ature |12| , po| , [21] , 22 1, while for sources with temperature gradients some numerical 



checks of the approximate analytical approximations developed in [13, 14, y^] have 
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been performed in p^, 24]. In |23| the validity range of the analytical formulae given 



in was determined. A comparison with experimental data from 200 A GeV/c S+Pb 
collisions |25| in [^] revealed, however, that the required source parameters lie outside 
this range of validity. A quantitative comparison with data thus requires numerical 
studies. 

In this paper we present a numerical parameter study of the influence of flow and 
temperature gradients on the correlation radii. One goal of this investigation is to 

24]: based on a series of approximations. 



settle a question which was left open in 
the authors of [0] had found a common l/^/M± scaling law for all three HBT radii 
Rs, Ro, Ri in a Cartesian parametrisation of the correlator. This was desirable in 
view of the S+Pb data of ]|2^ where agreement with such a common scaling law was 
reported (although supported only by three data points in each Cartesian direction). 
The data of l2^, 



28, 



\ 30, 31] from S+S, S+Ag, S+Au, and Pb+Pb collisions, on 
the other hand, all show a much stronger Mj^-dependence of the longitudinal HBT radii 
compared to the transverse ones. This is in qualitative agreement with the theoretical 
studies in ]|l^, 20, 21] which, for sources with constant temperature and boost-invariant 
longitudinal, but weaker transverse expansion, show a weaker M_|_-dependence for Rg 
than for Ri. In |12, 21] it was argued that, if one fits the M_|_-dependence of the HBT 
radii with a negative power law, Ri{M_[_) ~ M^°'' ]^ 

be proportional to the rate of expansion in direction i; this feature cannot be obtained 

We show 



^, [31[], the power a, itself should 
this feature cannot 

within the simple saddle point approximation used in ]10, 11, 13, 14, 16,|32 
here that, within the phenomenologically allowed parameter range, this remains true 
even in the presence of temperature gradients. 

In contrast to the work of ]^, 24] we also compute the HBT radii in the Yano- 
Koonin-Podgoretskii (YKP) parametrisation of the correlation function p3| , 34, 32] 
which have a more straightforward and simpler interpretation in space-time 
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than the HBT parameters from the Cartesian (Pratt-Bertsch-Chapman) parametrisa- 
tion ]10, As far as we know, the effects of temperature gradients on the YKP 
parameters have not been studied before. 

We want to stress that resonance decays are not addressed in this paper. It is known 
that they can strongly affect the correlation function ]35, 36, 37, 



3|, |g] in which 

case a Gaussian parametrisation as employed here becomes questionable ]^, 42]. 
Thus, while our calculations refer both to pion and kaon correlations, only the kaon 
results and those for high-K_L pions (where resonance decays can be neglected) can be 
directly related to data. For low-K± pion pairs our results show the features of the 
contribution from directly emitted pions to the correlator. 



2 Formalism 

For chaotic sources the two-particle correlation function is in very good approximation 
given by the formula J43, p5|] 



C{q,K) ~ 1 + 



]/(iS5(x,K)e*«-^] 
\J d'^xS{x,K)\'^ 



(1) 



Here S{x,K) is the emission function (single-particle Wigner phase-space density) of 
the source [|^, 44, 45], and K = ^{pi + P2) (with pi^2 on-shell) is the average pair 
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momentum while q = Pi — P2 denotes the momentum difference between the two 
particles. 

One usually parameterizes the correlation function by a Gaussian in q |20, 21 1. 
Since q satisfies the on-shell constraint q ■ K = 0, only three of its four components 
are independent. This leaves room for various (mathematically equivalent) Gaussian 
parametrisations, using different sets of independent g-components. Here we will con- 
sider only azimuthally symmetric sources and evaluate the parameters (HBT radii) of 
the Cartesian and of the Yano-Koonin-Podgoretskii (YKP) parametrisations, in the 
commonly used coordinate system where the z axis defines the beam direction and K 
lies in the x — z plane. The x axis is customarily labeled as o (for outward), y as s (for 
sideward), and z as I (for longitudinal). 

The Cartesian parametrisation of the correlator employs the three spatial compo- 



nents of q in the form [10, 11] 



C(q, K) = 1 + exp -RliK)ql - RliK)ql - Rf{K)qf - 2Rl{K)qoqi 



(2) 



The four K-dependent parameters (HBT radii) are then given by linear combinations 
of the space-time variances of the effective source of particles with momentum|^ K 

m, MM- 

r.2\ 



i?^(K) 
Rl{K) 

Rim 



{{x 
{{x 



-P±i)^), 
-I3±i){z 



■Pit)) 



(3) 
(4) 
(5) 
(6) 



The angular brackets denote space-time averages taken with the emission function 
S(x, K) of the effective source at momentum K, and coordinates with a tilde are 



measured relative to the center of the effective source 
In the LCMS (Longitudinally CoMoving System 

^ol 



R^i simplify to 



where Pi = Ki 



Rf 

rIi 



{{x - (5 J) ~z) 



-0, Rf and 

(7) 
(8) 



For further discussions of these expressions see |To| , 11 ] . 

The YKP parametrisation employs the components q^, qj_ = \/qi + qj, and qi in 
the form H, |^, |l| 



C(q,K) = 1 + exp 



-Rlmql 



^S(K) {ql 



0\2 



RliK)+Rl{K))iq.U{K)) 



(9) 

The YKP radii R±, and Rq are invariant under longitudinal boosts of the mea- 
surement frame; they measure (in some approximation) the transverse, longitudinal, 
and temporal extension of the effective source in its longitudinal rest frame (called 
Yano-Koonin (YK) frame) ||3^, ^]. This frame moves with the YK velocity v(K), 
defined by the fourth fit parameter U(K) via 

U{K) = 7(K) (1, 0, 0, ^;(K)), 7(K) = \ . (10) 

VI - ^^(K) 



^We denote the dependence on the pair momentum alternatively by K or if where in the latter 
case the on-shell approximation K'' = \/m?- ■\- is implied. 
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The YKP parameters can be calculated from the emission function through relations 



similar to (|)-(|): Defining |g, || 

A = {\t 

B = 

C = 

with = X + iy such that (^^) = {x^ — y^), they are given by 




A + B 



1 




2C 
A + B 



(11) 
(12) 

(13) 

(14) 

(15) 
(16) 
(17) 



Note that the Yano-Koonin velocity and thus R\\ and i?0) fire well-defined only for 
effective sources with {A + B)^ — 4C^ > 0. It will be shown in |1£] that this condition 



can be violated in particular for opaque sources. In the limit C — > the YK velocity 
V vanishes such that this condition defines the YK frame. In this frame R?, = B and 



A. 



As a side remark, we would like to note that the quantities A, B, C from Eqs. (11)- 
(|l^ ) are identical with the radius parameters introduced in 42], where a different 

Cartesian parametrisation based on the same g-components as in the YKP case was 
used: 



C(q,K) = 1 + exp[-i?i(K)gi - Rl{K)q^ - Rj{K){q^f - 2 Rlt{K)qiq^] 
One easily shows 



(18) 



Ri — A , i?, — B , R, 



-C . 



H — ^, ±1.^ — ±J , J-L^-i- — i_y . (19) 

Clearly, these parameters are not boost-invariant and thus don't lead directly to a 
simple space-time interpretation of the correlator. Their usefulness is rather of tech- 
nical nature in the actual data fitting procedure [p2| , ^] and, of course, the YKP 
parameters can be reconstructed from them via Eqs. (p!^)-([T^). 

Our calculations in the following sections of the correlation radii will be based on 
the expressions (D-d) and (0)-(0). 

3 Hydrodynamic parametrisation of the source 

For our studies we use slightly modified (see Appendix]^) model of p^ : 



Six,K)d\ 



M_|_ cosh(r/ — Y) 
(2^0^ 



exp 



K ■ u{x) 



T x) 



exp 



r'^ _ iv - Vo) 
'2i?2 2(Ar/)2 



^Please note that the first expressions given on the r.h.s of Eqs. (19b, c) in are only vahd if 
A + B>0. 
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r dr 

exp 



|2 



2(At) 



d-nrdrdcj). (20) 



This "emission function" parametrizes the distribution of points of last interaction in 
the source. The parametrisation (20) is motivated by hydrodynamical models with 



approximately boost-invariant longitudinal dynamics. It uses thermodynamic and hy- 
drodynamic parameters and appropriate coordinates; for a detailed discussion see, e.g., 
In the calculations we use the on-shell approximation Kq ~ Ek = \/vn? + as 
discussed in [Tl|] . 

The velocity field u{x) determines the dynamics of the source at freeze-out. We 



parametrize it by |12 



vF'ix) = (cosh 7/ cosh r/t(r), cos(/> sinh?7((r), sine/) sinh r/t(r), sinhr7cosh7y((r)) , (21) 

thereby implementing a boost-invariant longitudinal flow profile = z/t, with a linear 
radial profile of strength r\j for the transverse fiow rapidity^: 

m{r)=r]s^-. (22) 

Parameterizing K in the usual way through rapidity Y ^ transverse mass M_l and trans- 
verse momentum , the exponent of the Boltzmann term reads 

K ■ u{x) = M_L cosh(r/ — Y) cosh rit(r^ — sinh r7i(r) cos (/> . (23) 

For rjf = the emission function depends only on M± and not on Kj_. This scaling is 
broken by the transverse fiow. 

A special feature of the model suggested in IQ is the parametrisation of the tem- 
perature profile: 



r(x) To \^ 2Ry \ 2r2 

It introduces transverse and temporal temperature gradients which are scaled by the 
parameters a and d. Such a profile concentrates the production of particles with large 



M_|_ near the symmetry axis and close to the average freeze-out time [14, |I5|. Note 
that the space-time dependence of the temperature does not break the M_|_-scaling of 
the emission function in the absence of transverse fiow. 

The emission functions (effective sources for pions with given momenta) for different 
model parameters are shown as density contour plots in Figs. || and |2[ Transverse cuts 
of the emission function are displayed in Fig. 

Transverse flow is seen to decrease the effective source more in the sideward than 
in the outward direction (see Fig. ^3 and also Fig. 3 in |l^). Transverse temperature 
gradients, on the other hand, just reduce the homogeneity lengths without changing 
the shape of the source (cf. Figs. |l^,c). 

A temporal temperature gradient has no impact on the transverse source profile; 
its effect on the longitudinal profile is seen in Fig. |2| (top vs. bottom row) where it 

^The nonrelativistic approximation of this transverse profile used in [ pi] , p^ , p2[ is advantageous for 
analytic manipulations but not necessary if the correlation radii are evaluated numerically. 
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Figure 1: Transverse cuts of the emission function for midrapidity pions with transverse 
momentum K± = 500 MeV/c (in x-direction) , for the model parameters given in 
Table [l|. Left panel (a): neither flow nor temperature gradients; middle panel (b): 
transverse flow only with r]f = 0.7; right panel (c): transverse temperature gradient 
only with a = 0.8 



decreases the temporal width of the effective sources without otherwise changing the 
shape of the emission function. 

Figure ^ also shows the different shapes of the effective source in the center of mass 
system (CMS) of the fireball for midrapidity and forward rapidity pions. An interesting 
feature are the emission time distributions: one clearly sees that the boost-invariant 
longitudinal expansion and proper-time freeze-out leads to freeze-out at different times 
in different places and thus generates time distributions in a fixed coordinate system 
which cover a much larger region [12, 20 1 than the one corresponding to the intrinsic 
eigentime width encoded in the emission function (pO[). The strongly non-Gaussian 
shape of the longitudinal source distribution is known (T^ j to lead to non-Gaussian 
effects on the longitudinal correlator; as shown in |40| its Gaussian shape is, however, to 
some extent restored by resonance decay contributions (not considered here) which fill 
in the central region above the "hyperbola" in the thermal source shown in Figs. 

For our numerical model study we use, if not stated otherwise, the source (|2C 
the model parameters listed in Table ||. 




4 Influence of temperature and flow gradients on the 
HBT radii of a transparent source 

In this section we discuss the effects of various types of gradients in the source on the 
correlation radii. To be able to recognize their specific signals we first investigate the 
correlation radii in the absence of transverse flow and temperature gradients. At the 
end of this section we try with the help of temperature gradients to reproduce the 
l/\/M± scaling proposed in [0] for Rs and the emission duration. 

For the sake of clarity, let us list here the various reference frames which will appear 
in the following discussion. They differ by their longitudinal velocities. 

*We thank to Urs Wiedemann for making this point. 
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Figure 2: Contour plots of the z — t profiles of the emission function for pions. The 
profiles are calculated in the CMS. (a,b): pair rapidity y = 0, transverse pair momen- 
tum K^ = \ MeV/c. (c,d): F = 2, i^_L = 500 MeV/c. The sensitivity of the profiles 
to is weak. Source parameters as given in Table r/j = a = 0. Upper row (a,c): 
no temporal temperature gradient, d = 0; lower row (b,d): d = 1.5 



CMS Center of Mass System of the fireball. In this frame % = 0. 

LCMS Longitudinally Co-Moving System - a frame moving longitudinally with the 
particle pair, i.e., F = 0, /?; = 0. 

YK The Yano-Koonin frame. It moves longitudinally with the Yano-Koonin velocity 
V. The YKP radii measure the homogeneity lengths of the source in this frame. 

LSPS Longitudinal Saddle-Point System. This frame moves longitudinally with the 
point where the emission function for particles with a given momentum K has 
its maximum (point of maximal emissivity). 



Table 1: Values of model parameters used in numerical calculations 



temperature Tq 


100 MeV 


average freeze-out proper time tq 


7.8 fm/c 


mean proper emission duration Ar 


2 fm/c 


geometric (Gaussian) transverse radius K 


7 fm 


Gaussian width of the space-time rapidity profile Ar/ 


1.3 


pion mass 


139 MeV/c^ 


kaon mass vnj^^ 


493 MeV/c2 



7 



It is obvious that each particle pair (via its momentum K) defines its own LCMS, YK, 
and LSPS frame. The rapidity of the pair in the CMS will be denoted by Y^^ = Y — rjo- 

4.1 No temperature and flow gradients 

In Fig. ^we show the M± dependence of the correlation radii for both the Cartesian and 




200 400 600 800 1000 200 400 600 800 1000 



M^(MeV/c) M^(MeV/c) 

Figure 3: Correlation radii in the absence of temperature gradients and transverse 
flow, calculated for a transparent source with model parameters from Table ||. In the 
left (right) column the M± dependences of the YKP (Cartesian) radii are shown. The 
space-time variances from which they are calculated are evaluated in the LCMS. The 
HBT radii for pions (kaons) are shown as solid (dashed) lines, in three rapidity regions: 
^CM = 0> ^CM = 1-^ ^-iid Y^^ = 3 



Yano-Koonin-Podgoretskii parametrisations, for a source without transverse flow and 
temperature gradients. Non-trivial M± dependencies thus results solely from the effects 
of longitudinal expansion or have a kinematic origin. All HBT radii are calculated from 
space-time variances of the emission function evaluated in the LCMS. These curves are 
shown for later reference only; for a detailed interpretation of their features we refer 
to the existing literature |l2|, |3^, 20, 11]. Here we concentrate on a few relevant 
features: 

Due to the absence of transverse gradients, Rg 



R± is independent of K; it is not 
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affected by longitudinal gradients [^]. The outward Cartesian radius parameter Ro 
reflects the effective lifetime of the source (in the frame in which the expressions on 
the r.h.s. are evaluated) via [12, 47] 

Rl^ ^Rl- Rl = (3l(p) + 2 (xi) + lx^-f) . (25) 

Due to the absence of transverse flow the last two terms vanish here. The effective 
lifetime in the YK frame (i.e. in the rest frame of the effective source) is measured by 



Rq in the YKP parametrisation |2C, 21]. Due to different kinematic factors the lifetime 
effects are stronger and more easily visible in Rq than in R^. Figure Q also confirms 
the observation j^, ^] that the M±_ dependence of the effective emission duration is 
connected with the M±_ dependence of the longitudinal region of homogeneity which 
results from longitudinal expansion ]^]. 

The decrease of R\\ near the edge of the rapidity distribution results from an inter- 
play between the Gaussian space-time rapidity distribution and the Boltzmann factor 
in the emission function. It gives rise to a narrowing rapidity width of the effective 
source with increasing pair rapidity Y^^. At large M_|_, the longitudinal radii Ri and 
are independent of the pair rapidity Y^^^ , due to the longitudinal boost-invariance of 
the velocity profile in the Boltzmann factor which dominates the shape of the emission 
function in the limit M±_ oo. 

It is interesting to compare the longitudinal Cartesian radius parameter Ri in the 
LCMS frame with the longitudinal YKP radius For small M±, Ri is larger than 

although not by much, while at large M± the two parameters agree. This can be 
understood by recalling the relation 



Rf = (1 - Pif Rl + 7' (A - vf {rI + Rfi) (26) 



which, together with two other such relations ]20], expresses the mathematical equiv- 
alence of the Cartesian and YKP parametrisations. For large M± the source velocity 
V coincides with the longitudinal velocity (3i of the pair (which is zero in the LCMS). 
Then the second term in ( |2^ ) vanishes and Rf = R?, in the LCMS. For smaller values of 



M_|_ and Y^^ 7^ the pair and source velocities /5; and v are slightly different ]|2^; the 
resulting positive contribution from the second term in ( p6|) renders Rf > ii| unless 
Rq turns negative (cf. ]jl^). This feature was already observed by Podgoretskii ]34] 



who introduced the Yano-Koonin frame as the frame in which the production process 
is reflection symmetric with respect to the longitudinal direction. Within a class of 
non-expanding models he showed that in this frame the longitudinal source radius is 
minimal. On first sight this appears to contradict the laws of special relativity from 
which one might expect that the longitudinal homogeneity length of the source should 
be largest in the source rest frame and appear Lorentz contracted in any other frame. 
This argument neglects, however, the fact that different points of the homogeneity 
region generally freeze out at different times]^ (see Figs. ^,d). 

The cross term Rf^i shown in the lower right panel of Fig. |3| is required by the 
invariance of the correlation function under longitudinal boosts of the measurement 



frame [QG, 11]. Its value and M± dependence depends strongly on the measurement 



^We thank Jan Pisiit for a clarifying discussion on this point. 
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Table 2: The exponents found by a fit of the HBT radii with the function Ri{M±) oc 
i = l,\\,0, for a transparent source without transverse expansion and tempera- 
ture gradients 





Y 

^ CM 


1.5 3 


pions 
Jcaons 


0.568 0.585 0.626 
0.546 0.553 0.573 


pions 
" kaons 


0.568 0.542 0.478 
0.546 0.538 0.515 


pions 
kaons 


0.403 0.380 0.321 
0.202 0.197 0.185 



frame; for example, in the CMS the cross term is negative with smaher absolute value 
fi^ although its generic M± dependence remains similar. 

Except for Rq and R^i, all radius parameters shown in Fig. ^ scale with the trans- 
verse mass M± and show no explicit dependence on the particle rest mass. The curves 
for pions and kaons thus coincide. This scaling is broken for Rq and R^i by the ap- 
pearance of explicit factors P± in Eqs. @) and (^). A quantitative measure for the 
strength of the Mj_ dependence can be obtained by fitting the radius to a power law, 



Ri{M±) oc M_|_"% i = 1, 11,0 [27, The corresponding exponents are listed in Table 



|2[ The longitudinal radius parameters Ri and -Ry scale approximately with l/\/M± as 
predicted in [^. The scaling power is generally closer to -0.5 for kaons than for pions; 
this is due to their larger rest mass which reflects in larger values for M± where the 
saddle point approximation and the assumption of longitudinal boost-invariance in Q 
become better. For the temporal YKP parameter Rq the l/^/M± scaling law does not 
hold - it decreases more slowly, especially at larger M_|_ and for kaons. 

4.2 Effects from individual types of gradients 

In this subsection we investigate separately the effects of specific types of gradients in 
the emission function on the two-particle correlations, discussing at the same time the 
corresponding single particle spectra. A simultaneous analysis of one- and two-particle 
spectra is required for a clear separation of thermal and collective features of the source 
and for a complete reconstruction of its emission function |14, pO, 51]. 



We begin with the discussion of gradient effects on the single particle spectra. 
Figure ^ shows the influence of transverse and temporal temperature gradients and of 
transverse flow on the m_L-spectra at flxed Y,-,^ = 0. Qualitatively very similar features 
are seen at other rapidity values and in the rapidity integrated spectra. One sees that 
even very strong temporal and transverse temperature gradients do not cause any major 
effects on the single particle spectra, compared to the case without any gradients. Their 
main effect is a change of the normalisation because the colder regions of the emission 
function contribute less; their contribution is also concentrated at smaller momenta, 
leading to a somewhat steeper slope of the transverse mass spectra at low m±. 

This situation changes dramatically when the source develops transverse collective 
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Figure 4: Transverse mass spectra at midrapidity from a transparent source. Black 
lines correspond to pions, gray lines to kaons. In the left plot the effects from different 
types of gradients are shown individually, in the right plot we study the interplay of 
transverse temperature gradients with transverse flow 



flow. In the absence of temperature gradients, transverse flow simply flattens the m±- 
spectra |52]. For large ^ ttt-q (i.e. relativistic momenta) this can be described 

in terms of an effective blueshifted temperature T^g = T ^J^^^^ which is the same 

for all particle species while at small m± the flattening is even stronger and depends 
on the particle mass via the nonrelativistic relation Tefj = T + m{v{r))'^ ||5^. These 
features have now been clearly observed in the heavy collision systems studied at the 
Brookhaven AGS and CERN SPS [^]. Transverse flow clearly also breaks the m± 
scaling between pions and kaons: the pion and kaon spectra in the left panel of Fig. ^ 
have different slopes and normalisations. 

Additional temperature gradients which are imposed on top of the transverse ex- 
pansion flow affect the spectra as shown in the right panel of Fig. ^ A temporal 
gradient of the source temperature has no qualitative effect on the shape of the m± 
spectra and only affects their normalisation. A transverse spatial gradient of the tem- 
perature, however, interferes seriously with the transverse flow by strongly reducing its 
effect on the slope of the m± spectra. For a strong transverse temperature gradient 
with a = 5 as shown in the Figure, a comparison with the left panel shows that the 
flow effects become nearly invisible. An analytical discussion of this behaviour is given 
in where it is shown to occur if the transverse homogeneity length generated by 
temperature gradient becomes much smaller than the one generated by the transverse 
velocity gradient. It follows from this discussion that it is not possible to generate 
a strong (additional) M^-dependence of the transverse HBT radius parameter from 
transverse temperature gradients (see below) without at the same time reducing the 
flow effects on the single particle m± spectra. 

We now proceed to a discussion of the HBT radius parameters. The effect of 
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temporal temperature gradients is studied in Fig. ^ As expected, they affect mostly 
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Figure 5: The effects of temporal temperature gradients on the temporal YKP pa- 
rameter Rq (upper panel) and the difference — R^ in the Cartesian parametrisation 
(lower panel). Solid (dashed) lines are for pions (kaons). The calculations are done for 
midrapidity pairs, Y^^ = 0, and for a = r/j = 



those HBT parameters which are sensitive to the effective lifetime of the source, namely 
the temporal YKP parameter Rq and the Cartesian difference R^ — R^. These are 
reduced by increasing d. The origin of this effect can be seen in Fig. ^ which shows that 
larger values of d concentrate the high temperature region (with large emissivity) to a 
narrower region in the z-t plane. The narrowing occurs predominantly in the temporal 
direction, but affects slightly also the longitudinal homogeneity length. Going from 
central rapidity to pairs at forward rapidities, all curves change in a similar way as 
shown in Fig. |^. 

Transverse temperature gradients leave their traces only in the transverse homo- 
geneity length measured by Rg = R± ■ In Fig. ^ this is shown for pions and kaons at 
mid-rapidity. (Note that in our model Rg is rapidity independent.) 

The most important feature of both types of temperature gradients is that they do 
not break the M± scaling of the YKP correlation radii found in previous subsection. 
Note also that the temperature gradients have no qualitative impact on the YK rapidity 
y^j^ associated with the YK velocity v. 

Many of these features change in the presence of transverse flow. This is shown 
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Figure 6: The effect of a transverse temperature gradient: The transverse radius pa- 
rameter Rg = as a function of the transverse mass of the pair, for two different 
values of a. The coincidence of the solid (dashed) hues for pions (kaons) reflects the 
M_L-scaling of the transverse radius in this case, d = rjf = 0, Yf-.y^ = 



in Fig. ^ for pairs with rapidity Y^^^ = 1.5 in the CMS. All radius parameters have 
been calculated from space-time variances evaluated in the LCMS. The choice of a non- 
zero CMS pair rapidity ensures a non- vanishing cross term Roi- The generic rapidity 



dependence of the radii follows the behaviour discussed in Sec. 4.1. 

Since transverse flow introduces velocity gradients in the transverse direction, Rs{= 
R±) and Ro decrease with increasing rjf. Less obvious is, however, the rather strong 
effect of transverse flow on the "temporal" YKP parameter i?g. For higher transverse 
masses the latter even begins to increase with M±. This does not, however, reflect 
the behaviour of the emission duration which continues to decrease for particles with 
higher M±. The increase of Rq is rather caused by correction terms expressing the 
difference between R^ and the lifetime po| , 

Rl-{P) = ^{xi) + ^{x'-f). (27) 

Especially the second term on the r.h.s. grows appreciably with increasing r]f, and 
even more so for pions than for kaons. This was studied in detail in (cf. Fig. 2 of 
that paper) where a more detailed explanation of the curves for Rq shown here can be 
found. 

The most important feature of Fig. ^ is the breaking of the M± scaling of the YKP 
radii by the transverse flow. 

4.3 Interplay of all gradients 



In 1 14] it was claimed that an interplay of the various types of gradients discussed 



above can lead to a common l/\/M± scaling law for all three Cartesian HBT radii, Ri, 
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Figure 7: The correlation radii for a transparent source of constant temperature fea- 
turing transverse collective flow. The rapidity of the pairs was chosen as = 1.5. 
Left column: YKP radius parameters. Right column: Cartesian radius parameters. 
Solid (dashed) lines correspond to pions (kaons) 



Rs, and Ro, as well as for the effective lifetime of the source. The condition for such 
a common scaling is that the geometric extension is large compared to the "thermal 
length scales" generated by temperature and flow gradients |14|. The claim of |14| 
was based on approximate analytical expressions for these radius parameters which 
were obtained by evaluating the corresponding integrals over the emission function 
within an improved saddle-point approximation scheme. In this subsection we present 
a numerical test of this claim. To this end we select a set of parameters (a = 0.6, 
d = 10, rjf = 0.2) for which at mid-rapidity the conditions for such a scaling, as 
given in Q], are satisfied as much as possible, without leaving the phenomenologically 
realistic range. It should be mentioned, that our model differs slightly from that of |T^ . 
The differences are explicitly treated in Appendix ^ where they are shown to be small, 
and our findings are not affected by them. 

The numerically evaluated correlation radii for this parameter set are plotted with 
thick black lines in Fig. The analytical approximation according to |14] is shown 
by the thin lines for comparison. With thick grey lines the numerically computed 
results from original model of |14| introduced in Appendix ^ are plotted. We agree 
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Figure 8: Correlation radii in the Cartesian parametrisation for a transparent source 
with parameters a = 0.6, d = 10, ijf = 0.2. The thick black hues represent the exact 
result from a numerical evaluation of the space-time variances, the thick grey lines the 
result of numerical evaluation with the model of [0] with corresponding parameters 
(see Appendix |A|), the thin lines use the approximate analytical formulae given in |14|. 
Solid (dashed) lines are for pions (kaons). In cases of Rg and Ro thick black and grey 
lines overlap 



with the conclusions drawn from a similar comparison presented in |23| that for kaons 
the approximate analytical formulae agree with the numerical curves within 15%. For 
pions with transverse momenta below 600 MeV the discrepancies are larger because 
one of the conditions of validity of the analytical formulae (sufficiently large M±) 
is violated. However, the analytical approximation does not give at all the (slight) 
breaking of the M± scaling of Ri by the transverse flow, and at mid-rapidity it misses 
the initial rise of Ro at small Mj_. This last point in particular is a serious problem 
if one wishes to extract an estimate of the effective lifetime according to Eq. (^). 
Similar (dis-) agreement at a comparative level is seen for forward rapidity pairs at 
Y^^ = 1.5, with the exception of the cross term R^i which is strongly overestimated by 
the approximation. The analytical approximation of [0] is restricted by the condition 
that the flow rapidity of the point of maximum emissivity in the LCMS is small. Since 
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Table 3: The exponents ctj found in a fit of the numericahy computed correlation radii 
to a power law Ri oc Mj^"' , for a source with the parameters given in Fig. |^ 





Y 

^ CM 


1.5 3 


pions 

as 1 

kaons 


0.257 0.253 0.240 
0.306 0.304 0.297 


pions 

Oil , 

Jcaons 


0.512 0.527 0.565 
0.519 0.524 0.538 


pions 
" kaons 


0.512 0.496 0.453 
0.519 0.512 0.494 


pions 
kaons 


0.373 0.360 0.324 
0.225 0.223 0.262 



the latter coincides with good approximation with the YK rapidity [^] and, in our 
model, the YK rapidity is found to converge to Y in the limit M± — > oo |20, 21|, the 
point of maximum emissivity moves to zero LCMS rapidity in that limit. This explains 
why for large M±_ the analytical and numerical results agree. 

On a superficial level the analytical approximations of |14] thus don't seem to be 
doing too badly. We want to check, however, whether the results also confirm the 
suggested 1 /yj M^-scaling. In |^l| it was shown that in general the strength of the 
M_|_-dependence of Rg and Ri (resp. R^ and R\\) is correlated with the strength of 
collective flow in the transverse resp. longitudinal directions. One way of quantifying 
the strength of the M^-dependence of the correlation radii is to fit them to a power 
law Ri (X Mj^"' [i = s,l, ||,0) and to study the values of the exponents a^. According 
to the approximate formulae of |jl^, at mid-rapidity all these powers should be equal 
to 0.5. The Qfj-values obtained from the numerically evaluated correlation radii are 
listed in Table |3|. One sees that only the longitudinal radius parameters (both in YKP 
and Cartesian parametrisations) follow an approximate M± scaling law, similar to 
previous studies without temperature gradients [|l2|, It reflects the boost-invariance 
of the longitudinal expansion velocity profile |^]. The transverse and temporal radius 
parameters Rg = R± and Rq scale much more weakly with M±. For mid-rapidity 
kaons Rq provides a good estimate for the emission duration because for T]f = 0.2 the 
correction terms on the r.h.s. of (27) are small; thus we find that also the effective 
lifetime of the source does not scale with M± . Thus, while the conditions for the 
saddle point approximation in |14] are satisfied with sufficient accuracy by our (semi- 
realistic) parameter set, the stronger conditions required for the common 1/^/ M±- 
scaling are not. 



5 Conclusions 

Let us shortly summarize the most important results: 

Transverse temperature gradients and transverse flow have similar effects on the 
transverse HBT radii Rs = R±, but not on the single particle spectra. Both lead to a 
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decrease of Rs with increasing M± , but while transverse flow flattens the single particle 
spectra, a transverse temperature gradient only reduces their normahsation. A weak 
M_L-dependence of Rg due to a moderate transverse flow can be increased by adding 
a transverse temperature gradient, but only at the expense of simultaneously reducing 
the flow effects on the single particle spectra. This is important for the interpretation 
of experimental data. 

For transparent sources, transverse temperature gradients preserve the M_|_-scaling 
of the YKP radius parameters while transverse flow breaks it. This breaking of Mj_- 
scaling is weak, however, the most sensitive parameter being Rq. 

Temporal temperature gradients affect mostly the temporal HBT radii, i.e. the 
difference R^—R^ in the Cartesian parametrisation and Rq in the YKP parametrisation. 

We could not confirm the existence of a common l/\/M± scaling law for all three 
Cartesian HBT radius parameters and the effective source lifetime in the case of the 
"thermal lengths" being smaller than the geometric lengths. For realistic model pa- 
rameters, the M_|_-dependence for Rg and Rq is always much weaker than for the 
longitudinal radius parameters Ri resp. Rii. 
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A The model of Csorgo and Lorstad 



In this Appendix we point out the differences between our model and that of [14| and 
study their impact on our findings from Sec. 4.3. 

The model of |14] is given by the emission function 



S{x,K)d^x 



■J^^M± cosh(r/ - Y) exp 



K ■ u{x) ^ /i(x) 



; , ,v exp 



2 At' 



T{x) T{x) 
dr] r d(l)dr . 



(28) 



Since in the used coordinates the measure reads d^x = dr t dr] r dcj) dr , it is seen that 
here another geometric eigentime distribution is used, namely 



H(r) = ^ 



1 



T V2tt At'^ 



exp 



2 At" 



(29) 



This differs from our distribution by the pre- factor Tq/t. For the comparison with our 
model, the values of Tq and Ar' are obtained via the saddle point approximation to 
our Gaussian eigentime distribution multiplied by r from the Jacobian. This gives 



■^0 
Ar' 



-^0 + \lr?: + 4At2 



r^Ar 



At2 + 



(30) 
(31) 
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The remaining geometry of the model ( pq ) is encoded in chemical potential 

T{x) To 2 m 2A??2 ' ^"^ > 

so this is the same as in our model (pO|). There is a slight difference in the formulation 
of temperature gradients: 



rr., . ^i-'- o^^^+d- (33) 

r(x) To \ 2t'V \ 2r'2 ^ ' ^ ^ 



where the transverse gradient scales with c'/tq (and not with a/R). Furthermore, the 
transverse flow rapidity is given as 



rit = ri'f^, (34) 

instead of rjfr/R. This requires appropriate re-scaling of a' and tj'jt. Also, d' has to 
be re-scaled relatively to d by the factor Tq/tq. Note finally that in the numerical 
calculation we have used the expansion four-velocity field as given by the relativistic 
formula ( [2l| ) with the transverse rapidity profile (p4[), while in [0] the non-relativistic 
approximation to the transverse expansion is used in the formulation of the model. 

The results are plotted with grey lines in Fig. |8|. It is clearly seen that the difference 
between the mentioned two models is small. 

For the computation of the curves using the analytical formulae found in (thin 
lines in Fig. ^ we used the re-scaled model parameters as explained above. 

To be in the region where the l/V M_L-scaling should take place we have chosen: 
a = 0.6, d = 10, rjf = 0.2. Note that with the appropriately re-scaled parameter values 



(see above) we satisfy the conditions given in |Q which need to be satisfied for a 15- 
20% agreement between the numerical calculation and the analytical approximations 
of the model (|2J 



To be sure that the slope parameters found in Table are not an artefact of the 



difference between our model and that of |14] we did the same fit with the numerically 



calculated radii resulting from (^8|) and even with the analytically computed curves. 
Results of that fit are listed in Table ^. The differences of the exponents listed in 
Tables | and | are in case of as, ai and ay within 1%, ao's differ up to 10%. The 
qualitative discussion concerning the exponents from Sec. |4.3| remains valid for the set 
given in Table The same is true for the analytically determined radii. 

This again supports the conclusion, that even if good agreement between numerical 
and analytical results is achieved, for a realistic set of parameters the conditions for 

— 1/2 

common scaling of all three Cartesian radii with M , are not fulfilled. 



References 

[1] R. Hanbury Brown and R.Q. Twiss, Phil. Mag. 45, 633 (1954); 
and Nature 177, 27 (1956); 178, 1046 (1956); 178, 1447 (1956) 

[2] G. Goldhaber, S. Goldhaber, W. Lee, and A. Pais, Phys. Rev. 120, 300 (1960) 

[3] M. Gyulassy, S.K. Kauffmann, and L.W. Wilson, Phys. Rev. C20, 2267 (1979) 



18 



Table 4: The exponents ctj found in a fit of the numericahy computed correlation radii 
to a power law Ri oc MJ"' , for a source from 28 and with the parameters correspond- 
ingly re-scaled to approximate those used in Fig. ^. In brackets are the results of the 
same fit done with curves computed using the analytical approximation from |14] 





Y 

^ CM 


1.5 3 


pions 
as , 

kaons 


0.257 (0.301) 0.253 (0.274) 0.240 
0.306 (0.370) 0.304 (0.364) 0.297 


pions 
kaons 


0.511 (0.569) 0.526 (0.562) 0.564 
0.518 (0.531) 0.523 (0.535) 0.538 


pions 
" kaons 


0.511 0.495 0.452 
0.518 0.512 0.493 


pions 
kaons 


0.395 0.381 0.344 
0.245 0.242 0.279 



K 

[6: 
[7; 

[9 

[lo: 
[11 
[12: 

[13 
[14 

[15: 

[16 

[17: 



B. L5rstad, Int. J. Mod. Phys. A12, 2861 (1989) 

D. Boal, C.K. Gelbke and B. Jennings, Rev. Mod. Phys. 62, 553 (1990) 

W.A. Zajc, in: Particle Production in Highly Excited Matter, edited by 
H.H. Gutbrod and J. Rafelski, (Plenum Press, New York, 1993), p. 435 

S. Pratt, in: Quark-Gluon Plasma 2, edited by R.C. Hwa, (World Scientific Publ. 
Co., Singapore, 1995), p. 700 

U. Heinz, in: Correlations and Clustering Phenomena in Subatomic Physics, 
edited by M.N. Harakeh, O. Scholten, and J.H. Koch, NATO ASI Series B, 
(Plenum, New York, 1997) (Los Alamos eprint archive [nucl-th /9609029D 

A.N. Makhlin and Y.M. Sinyukov, Z. Phys. C39, 69 (1988) 

S. Chapman, P. Scotto, and U. Heinz, Phys. Rev. Lett. 74, 4400 (1995) 

S. Chapman, P. Scotto, and U. Heinz, Heavy Ion Physics 1, 1 (1995) 

U.A. Wiedemann, P. Scotto and U. Heinz, Phys. Rev. C53, 918 (1996) 

T. Csorgo and B. Lorstad, Nucl. Phys. A590, 465c (1995) 

T. Csorgo and B. Lorstad, Phys. Rev. C 54, 1390 (1996) 

T. Csorgo and B. Lorstad, in Proceedings of XXV International Symposium on 
Multiparticle Dynamics, Sept. 12.-16. 1995, Stara Lesna, Slovakia, edited by 
D. Bruncko, L. Sandor and J. Urban (World Scientific, Singapore, 1996), p. 661 

S.V. Akkelin and Y.M. Sinyukov, Z. Phys. C72, 501 (1996) 

H. Heiselberg and A. P. Vischer, Los Alamos eprint archive nucl-th /9609022 



19 



[18] H. Heiselberg and A. P. Vischer, Los Alamos eprint archive nucl-th/970303C 



[19] B. Tomasik and U. Heinz, continuation of this paper containing study of opaque 
sources - work in progress 

[20] U. Heinz, B. Tomasik, U.A. Wiedemann, and Y.-F. Wu, Phys. Lett. B382, 181 
(1996) 

[21] Y.-F. Wu, U. Heinz, B. Tomasik, and U.A. Wiedemann, Los Alamos eprint 



archive |nucl-th/9607044| , Z. Phys. C 77 (1997) in press 



[22] B. Tomasik, U. Heinz, U.A. Wiedemann, and Y.-F. Wu, Acta Phys. Slovaca 47, 
81 (1997) 

[23] T. Csorgo, P. Levai, and B. Lorstad, Acta Phys. Slovaca 46, 585 (1996) 

[24] T. Csorgo and B. Lorstad, Heavy Ion Physics 4, 221 (1996) 

[25] NA44 Coll., H. Beker et al., Phys. Rev. Lett. 74, 3340 (1995) 

[26] NA35 Coh., T. Alber et al., Z. Phys. C 66, 77 (1995); 

NA35 Coll., T. Alber et al., Phys. Rev. Lett. 74, 1303 (1995) 

[27] T. Alber, PhD thesis, MPI fiir Physik, Miinchen (1995), unpubhshed 

[28] T. Alber for the NA35 and NA49 Coh., Nucl. Phys. A 590, 453c (1995) 

[29] K. Kadija for the NA49 Coll., Nucl. Phys. A 610, 248c (1996) 

[30] S. Schonfelder, PhD thesis, MPI fiir Physik, Miinchen, 1996, unpublished 

[31] H. Appelshauser, PhD thesis, Johann Wolfgang Goethe-Universitat 
Frankfurt/Main, 1996, unpublished 

[32] S. Chapman, J.R. Nix, and U. Heinz, Phys. Rev. C52, 2694 (1995) 

[33] F. Yano and S. Koonin, Phys. Lett. B78, 556 (1978) 

[34] M.I. Podgoretskii, Sov. J. Nucl. Phys. 37, 272 (1983) 

[35] J. Bolz et al., Phys. Lett. B300, 404 (1993); and Phys. Rev D47, 3860 (1993) 

[36] B.R. Schlei et al., Phys. Lett. B376, 212 (1996) 

[37] U. Ornik et al., Los Alamos eprint archive |hep-ph / 9604323 

[38] H. Heiselberg, Phys. Lett. B379, 27 (1996) 

[39] T. Csorgo, B. L5rstad, and J. Zimanyi, Z. Phys. C71, 491 (1996) 



[40] U.A. Wiedemann and U. Heinz, Los Alamos eprint archive |nucl- th /9611031 
submitted to Phys. Rev. C 



[41] U.A. Wiedemann and U. Heinz, Los Alamos eprint archive |nucl-th/9610043 



[42] U.A. Wiedemann and U. Heinz, Acta Phys. Slovaca 47, 95 (1997) 

20 



[43] E. Shuryak, Phys. Lett. B44, 387 (1973); Sov. J. Nucl. Phys. 18, 667 (1974) 

[44] S. Pratt, Phys. Rev. Lett. 53, 1219 (1984); Phys. Rev. D33, 1314 (1986) 

[45] S. Chapman and U. Heinz, Phys. Lett. B340, 250 (1994) 

[46] M. Herrmann and G.F. Bertsch, Phys. Rev. C51, 328 (1995) 

[47] T. Csorgo and S. Pratt, in Proceedings of the Workshop on Relativistic Heavy 
Ion Physics at Present and Future Accelerators, Budapest, 1991, edited by T. 
Csorgo et al. (MTA KFKI Press, Budapest, 1991), p. 75 

[48] J.D. Bjorken, Phys. Rev. D27, 140 (1983) 

[49] U. Heinz, Nucl. Phys. A610, 264c (1996) 

[50] S. Chapman and J.R. Nix, Phys. Rev C54, 866 (1996) 

[51] U. Heinz, B. Tomasik, and U.A. Wiedemann, work in progress 

[52] K.S. Lee, U. Heinz, and E. Schnedermann, Z. Phys. C48, 525 (1990); 

E. Schnedermann, J. Sollfrank, and U. Heinz, Phys. Rev. C48, 2462 (1993) 

[53] Proceedings from Quark Matter '96, (Heidelberg, 20.-24.5.1996), edited by P. 
Braun-Munzinger et al., Nucl. Phys. A 610 (1996) 



21 



